Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
神川 豊; 鈴木 真琴; 安掛 寿紀; 村上 貴彦; 森田 祐介; 椎名 秀徳; 福島 学; 平根 伸彦; 大内 靖弘
JAEA-Technology 2023-030, 57 Pages, 2024/03
航空機落下事故に関するデータが原子力規制庁により更新されたことに伴い、原子力科学研究所における航空機落下確率を再評価するため、経済産業省原子力安全・保安院「実用発電用原子炉施設への航空機落下確率の評価基準について(内規)」に基づき評価を行い、原子力科学研究所の各施設における航空機落下確率を評価した。評価の結果、航空機落下確率の総和は最大となる放射性廃棄物処理場において 5.6810回/(炉・年))であり、航空機落下を「想定される外部人為事象」として設計上考慮する必要があるか否かの基準である 10回/(炉・年))を超えないことを確認した。
渡辺 勢也*; 河原 淳*; 青木 尊之*; 杉原 健太; 高瀬 慎介*; 森口 周二*; 橋本 博公*
Engineering Applications of Computational Fluid Mechanics, 17(1), p.2211143_1 - 2211143_23, 2023/00
被引用回数:1 パーセンタイル:61.52(Engineering, Multidisciplinary)津波氾濫や大雨の斜面災害では、多くの浮遊物や流木、丸太などが流れに含まれる。剛体衝突による構造物の被害は、水圧による被害よりも深刻である。浮遊物体を含む自由表面流れを研究するためには、大規模計算が可能で高性能な自由表面流のシミュレーションコードの開発が必要となる。本研究ではキュムラント格子ボルツマン法と粒子ベースの剛体シミュレーションを組み合わせた単相自由表面モデルを提案する。剛体間の接触相互作用は個別要素法で計算される。解析精度の向上と計算の高速化のために8分木ベースの局所細分化格子法を導入し、自由表面や固体表面近傍には高解像度の格子を割り当てた。提案モデルの精度検証のために、八戸工業大学と神戸大学の15m水槽と70m水槽で2種類の津波実験を行った。シミュレーションの結果、漂流速度、捕捉木片数、積層角について実験と良い一致を示した。
市原 義孝*; 中村 尚弘*; 森谷 寛*; 堀口 朋裕*; 崔 炳賢
日本原子力学会和文論文誌, 21(1), p.1 - 14, 2022/03
本研究は、鉄筋コンクリート構造物の非線形性の影響を近似的に等価線形解析手法による地震応答解析で評価することを目的に、1996年にOECD/NEAによる国際解析コンペで使用された原子炉建屋耐震壁終局応答試験の三次元有限要素法によるシミュレーション解析を実施した。耐震壁の等価剛性及び等価減衰は、日本電気協会が提案するトリリニア型スケルトンカーブ、Cheng et al.が提案する履歴曲線より求め、せん断ひずみ調整ファクターは感度解析より0.70に決定した。その結果、せん断ひずみ=2.010程度までの試験体上部の卓越振動数,最大応答加速度,最大応答変位,慣性力-変位関係,床応答スペクトルを良く再現できることを明らかにした。本報における等価線形解析は、=4.010程度の終局破壊時の最大応答変位を過小評価している。このため、破壊直前の急激な変位の増大を含む試験結果の評価に本手法を適用する場合は、その適用性に十分留意する必要がある。
長谷川 雄太; 今村 俊幸*; 伊奈 拓也; 小野寺 直幸; 朝比 祐一; 井戸村 泰宏
Proceedings of 13th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Heterogeneous Systems (ScalAH22) (Internet), p.10 - 17, 2022/00
格子ボルツマン法(LBM)に基づく数値流体力学シミュレーションおよび局所アンサンブル変換カルマンフィルタ(LETKF)によるアンサンブルデータ同化をNVIDIA A100 GPU搭載スパコンに対して実装し、および最適化した。LBMとLETKFの協働のため、データ転置通信を実装し、LETKFのデータ依存性に基づいて計算,ファイルI/O、および通信のオーバーラップにより通信を最適化した。2次元等方乱流,アンサンブル数,格子点数の条件において、通信を最適化した実装は、LETKFを並列化しない単純な実装に対して3.85倍の高速化を達成した。LETKFの主要な計算カーネルはの実対称密行列の固有値分解であり、その計算のため、バッチ形式固有値分解ソルバEigenGを新たに開発した。EigenGによるバッチ形式固有値分解は、既存ライブラリであるcuSolverに対して64倍の高速化を達成した。
長谷川 雄太; 青木 尊之*; 小林 宏充*; 井戸村 泰宏; 小野寺 直幸
Parallel Computing, 108, p.102851_1 - 102851_12, 2021/12
被引用回数:2 パーセンタイル:32.94(Computer Science, Theory & Methods)GPUスーパコンピュータに対して格子ボルツマン法(LBM: lattice Botltzmann method)およびforest-of-octreesに基づくブロック構造型の局所細分化格子(LMR: local mesh refinement)を用いた空力解析コードを実装し、その性能を評価した。性能評価の結果、従来の空間充填曲線(SFC; space-filling curve)に基づく領域分割アルゴリズムでは、本空力解析において袖領域通信のコストが過大となることがわかった。領域分割の改善手法として本稿では挿し木法を提案し、領域分割の局所性とトポロジーを改善し、従来のSFCに基づく手法に比べて通信コストを1/31/4に削減した。強スケーリング測定では、最大で1.82倍の高速化を示し、128GPUで2207MLUPS(mega-lattice update per second)の性能を達成した。弱スケーリング測定では、8128GPUで93.4%の並列化効率を示し、最大規模の128GPU計算では44.73億格子点を用いて9620MLUPSの性能を達成した。
広田 憲亮; 柴田 裕司; 武内 伴照; 大塚 紀彰; 土谷 邦彦
Journal of Nuclear Science and Technology, 57(12), p.1276 - 1286, 2020/12
被引用回数:0 パーセンタイル:0.01(Nuclear Science & Technology)ケーブル長に沿った電位分布の安定性を達成することを目的として、高温条件に供された時の電気的特性に対する無機絶縁(MI)ケーブルの材料の影響を伝送テストによって調査した。その結果、MIケーブルの絶縁材料として、酸化アルミニウム(AlO),酸化マグネシウム(MgO)が選定され、ケーブルに沿った電圧降下の発生を確認した。有限要素法(FEM)に基づいた解析を実行し、終端部で検出された電位の漏れを評価した。伝送テストと解析による電圧降下の収率は、AlOおよびMgO材料のMIケーブルにおいてよい一致を示し、FEM解析結果と実験結果との相対的な関係を再現した。電圧降下を抑えるため、同様のFEM解析を行い、芯線直径()と芯線間距離()を変化させた。の変化を考えた場合、MIケーブルの電位分布は、絶縁材料の直径()をで割って得られる比率が0.35で最小電圧降下となった。を変化させた場合、最小電圧降下は0.5のl/であった。
杉田 裕; 谷口 直樹; 牧野 仁史; 金丸 伸一郎*; 奥村 大成*
日本原子力学会和文論文誌, 19(3), p.121 - 135, 2020/09
使用済燃料を直接処分するための処分容器の一連の構造解析を実施して、処分容器の必要な耐圧厚さの予察的な検討結果を示した。直線,三角形,正方形に配置された2, 3, 4体の使用済燃料集合体を収容するように処分容器を設計した。処分容器の胴体部分および蓋部分の必要な耐圧厚さを評価するため、使用済燃料集合体の収容スペースの離間距離をパラメータとした。この検討では、応力評価ラインの設定の妥当性や解析におけるモデル長の影響など、解析に関する技術的知識も得られた。そして、これらは、さまざまな条件下で同様の評価を実行したり、より詳細な評価を進めたりするための基盤として参考となるものである。
戸栗 智仁*; 矢萩 良二*; 沖原 光信*; 竹内 伸光*; 黒崎 ひろみ*; 松井 裕哉
JAEA-Technology 2018-017, 161 Pages, 2019/03
日本原子力研究開発機構東濃地科学センターでは、「国立研究開発法人日本原子力研究開発機構の中長期目標を達成するための計画」(平成27年4月1日平成34年3月31日)に基づき、三つの必須の課題(地下坑道における工学的対策技術の開発、物質移動モデル化技術の開発、坑道埋め戻し技術の開発)についての調査研究を進めている。本報告書は、これらの必須の課題のうち、坑道埋め戻し技術の開発として計画していた瑞浪超深地層研究所研究坑道を利用した坑道一部埋め戻し試験の全体計画等を策定した結果を述べたものである。
吉田 啓之; 上澤 伸一郎; 堀口 直樹; 宮原 直哉; 小瀬 裕男*
Proceedings of 11th Korea-Japan Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS-11) (Internet), 6 Pages, 2018/11
The radioactive aerosol removal equipment is used as one of the safety systems of nuclear reactors. In this equipment, micro particles of aerosol remove through gas-liquid interfaces of two-phase flow. The mechanism related to the removal of micro particles through the gas-liquid interface is not clear, a numerical evaluation method of performance of aerosol removal equipment is not realized. Then, we have started to construct a numerical simulation method to simulate removal of micro particles through gas-liquid interfaces. In this simulation method, detailed two-phase flow simulation code TPFIT is used as the basis of this method. TPFIT adopts an advanced interface tracking method and can simulate interface movement and deformation directly. In addition, to simulate the movement of particles, the Lagrangian particle tracking method is incorporated. By combining the interface tracking method and the Lagrangian particle tracking method, the interaction between interfaces and micro particles can be simulated in detail. To solve the Lagrangian equations of particles, fluid properties and fluid velocity surrounding aerosol particles are evaluated by considering the relative position of particles and gas-liquid interface, to simulate particle movement near the interface. In this paper, outline and preliminary results of this simulation method are shown.
守田 幸路*; 小川 竜聖*; 時岡 大海*; Liu, X.*; Liu, W.*; 神山 健司
Proceedings of 12th International Topical Meeting on Nuclear Reactor Thermal-Hydraulics, Operation and Safety (NUTHOS-12) (USB Flash Drive), 11 Pages, 2018/10
EAGLE炉内ID1試験は日本原子力研究開発機構によって実施され、FAIDUSと称される内部ダクト付き燃料集合体からの早期燃料流出を模擬したものである。試験で生じた早期ダクト破損は、燃料とスティールから構成される溶融プールからの高い熱流束によるものと解釈されている。試験後の分析からは、壁面に燃料クラストが形成されない状況において、高い熱伝導度を有するプール中の溶融スティールによって溶融プールからダクトへの伝熱が効果的に促進されたことが示唆されている。本研究では、多成分多相流の粒子法に基づいた完全ラグランジェ法を用いて溶融プールからダクト壁への熱伝達機構を分析した。プール中の溶融スティールと燃料の混合と分離挙動およびこれらの挙動がプールからダクトへの伝熱に与える影響を調べるため、燃料ピンの崩壊、溶融プールの形成およびダクト壁の破損に至る一連の挙動を模擬した。現在の2次元粒子法シミュレーションでは、10MW/mを超える壁面への大きな熱負荷は、核発熱を伴う液体燃料が壁面へ直接接触することによるものであることが示された。
岩元 大樹; 西原 健司; 八木 貴宏*; Pyeon, C.-H.*
Journal of Nuclear Science and Technology, 54(4), p.432 - 443, 2017/04
被引用回数:20 パーセンタイル:88.58(Nuclear Science & Technology)To investigate the applicability of the pulsed neutron source (PNS) method using a pulsed spallation neutron source (PSNS) for an on-line subcriticality monitoring system for an accelerator-driven system (ADS), a subcriticality experiment is conducted using Kyoto University Criticality Assembly (KUCA) in combination with the fixed-field alternating gradient (FFAG) accelerator. Reactivity values obtained from different traditional techniques, the area-ratio method and the -fitting method, are discussed with respect to the applicability to on-line subcriticality monitoring. The result shows that the area-ratio method robustly and accurately monitors subcriticality in shallow subcritical states of negative reactivity of up to a few dollars; however with this method, it faces problems with temporal fluctuations, spatial dispersion, and sensitivity to the proton-beam current with increasing depth of subcriticality. As a complement to this method, it is shown that the -fitting method alleviates such problems in deep subcritical state. Moreover, a proposed fitting technique using the maximum-likelihood estimation method based on the Poisson distribution is robust enough to be applicable for measuring negative reactivity of up to roughly nine dollars.
堀口 直樹; 吉田 啓之; 阿部 豊*
Nuclear Engineering and Design, 310, p.580 - 586, 2016/12
被引用回数:9 パーセンタイル:64.88(Nuclear Science & Technology)So celled multi Venturi scrubbers systems are installed in light water reactor systems to realize filtered venting at accident. This system is composed of many Venturi scrubbers (VS) and a water pool. However, the scrubbing mechanism of the self-priming VS including effects of gas mass flow rate and shape of the VS are understood insufficiently in the previous studies. Therefore, we started numerical and experimental study to understand the detailed two-phase flow behavior in the VS. In this paper, to understand the VS operation characteristics for the filtered venting, we performed numerical simulations of two-phase flow behavior in the VS. In the first step of this study, we perform numerical simulations of supersonic flow by the TPFIT to validate the applicability of the TPFIT for high velocity flow like flow in the VS. In the second step, numerical simulation of two-phase flow behavior in the VS including self-priming phenomena. As the results, dispersed flow in the VS was reproduced in the numerical simulation, as same as the visualization experiments.
矢野 直峰*; 山田 太郎*; 細谷 孝明*; 大原 高志; 田中 伊知朗*; 日下 勝弘*
Scientific Reports (Internet), 6, p.36628_1 - 36628_9, 2016/12
被引用回数:11 パーセンタイル:61.6(Multidisciplinary Sciences)We developed and employed a profile fitting method for the peak integration of neutron time-of-flight diffraction data collected by the IBARAKI Biological Crystal Diffractometer (iBIX) at the Japan Proton Accelerator Research Complex (J-PARC) for protein ribonuclease A and alpha-thrombin single crystals. Gaussian convolved with two back-to-back exponentials was selected as the most suitable fitting function, and a profile fitting algorithm for the integration method was developed. It was clearly demonstrated that the profile fitting method provides more accurate integrated intensities and model structures than the summation integration method at higher resolution shells.
中島 憲宏; 西田 明美; 宮村 浩子; 飯垣 和彦; 沢 和弘
可視化情報学会誌(USB Flash Drive), 36(Suppl.2), 4 Pages, 2016/10
組立構造を意識した有限要素解析により、原子力プラント全体での俯瞰的な耐震裕度評価と各部ごとの詳細な評価技術を実現する手段として、FIESTA(Finite element analysis for structure of assembly)と呼ぶコードを開発するとともに、組立構造解析技術の研究を進めている。本報では、観測値と計算解の照合を具体事例で例証した結果の可視化技術について報告する。地震波としては、震源地でマグニチュード5以上の7波を使い、「京」コンピュータにより地震応答解析を実施した。結合部を有する複雑構造物である原子力プラントの振動解析/耐震性評価過程を、Visual Analytics手法により実施した内容について報告する。これにより「組立構造解析」で耐震性評価に不可欠な計算結果の「確かさ」を向上する方法論をVisual Analytics手法により提案できた。
中島 憲宏; 西田 明美; 川上 義明; 鈴木 喜雄; 沢 和弘; 飯垣 和彦
Proceedings of 23rd International Conference on Nuclear Engineering (ICONE-23) (DVD-ROM), 5 Pages, 2015/05
数値解析結果の妥当性を推定する手段の一つを提案する。有限要素法による構造解析をスーパーコンピュータの並列計算機能を活用して実行し、解析結果の相違を分析するとともに、入力データである有限要素分割の粗密を変更した解析結果を分析し、数値計算的な確かさを推定するシステムを提案する。解析結果を表現する解析結果モデルの形成過程は、データベースで機能IDとその計算手順のリストによって記述する。解析モデルマネージャは、計算手順を記述したリストの順序によって、すなわち複数の数値計算手続きにより、シミュレーションを実行することで、目的とするシミュレーションの数値計算解を出力する。その結果、目的とするシミュレーション結果が複数生成されることから、これらの結果の相違を比較し分析することで、解の正確さを推定する。本論での数値実験は静解析と動解析で実施し、その正確さを判断するための必要な手続きを明らかにした。数値実験は、京を用いて行った。
加藤 由幹; 吉田 啓之; 横山 諒太郎*; 金川 哲也*; 金子 暁子*; 文字 秀明*; 阿部 豊*
Proceedings of 23rd International Conference on Nuclear Engineering (ICONE-23) (DVD-ROM), 8 Pages, 2015/05
In this study, visualization experiment for the bubbly flow in a horizontal pipe excited by oscillation acceleration was performed to understand bubbly flow behavior under earthquake acceleration. Liquid pressure was also measured at upstream and downstream of the test section. In addition, to consider detailed effects of pressure gradient on bubble motion, numerical simulation of two-phase flow in horizontal pipe with vibration was performed by a detailed two-phase flow simulation code with an advanced interface tracking method: TPFIT. Based on observed images and calculated results, bubble velocity was evaluated. It was confirmed that the pressure gradient amplitude increased with the increase of the frequency of the table. In addition, it was observed that the bubble velocity amplitude also increases with the increase of the frequency of the table. It was concluded that the bubble motion was strongly affected by the pressure gradient in the test section.
吉田 啓之; 永武 拓; 高瀬 和之; 金子 暁子*; 文字 秀明*; 阿部 豊*
Mechanical Engineering Journal (Internet), 1(4), p.TEP0025_1 - TEP0025_11, 2014/08
In this study, to develop the predictive technology of two-phase flow dynamics under earthquake acceleration, a detailed two-phase flow simulation code with an advanced interface tracking method TPFIT was expanded. In addition, the bubbly flow in a horizontal pipe excited by oscillation acceleration and under the fluctuation of the liquid flow was simulated by using the modified TPFIT. In the results, it was confirmed that the modified TPFIT can predict time dependent velocity distribution around the bubbles and shapes of bubbles qualitatively. The main cause of bubble deformation observed is large shear stress at the lower part of the bubble, and this large shear stress is induced by the velocity difference between the liquid phase and bubble. Moreover, we discussed about the difference between both effects of flow rate fluctuation and structure vibration on two-phase flow. In the results, bubble acceleration of the structure vibration case was larger than that of the flow rate fluctuation case. Finally, it was concluded that unsteady shear stress induced by vibration of the pipe wall was one of the main driving forces of bubble motion in structure vibration case.
Wang, X.*; Hou, Z.*; 池田 隆司; 寺倉 清之*
Journal of Physical Chemistry C, 118(25), p.13929 - 13935, 2014/06
被引用回数:10 パーセンタイル:33.24(Chemistry, Physical)グラフェンの端と欠陥での可能な窒素含有部のNMR化学シフトを第一原理計算により調べた。我々の計算結果はピリジン様窒素とグラファイト様窒素がN NMRにより容易に同定できることを示しており、これは実験と一致している。一方、ピリジニウム様窒素とピロール様窒素を識別することはこれらのN核のNMRシグナルが重なるために困難である。しかしながら、我々はシミュレーションからH NMRがこれらを区別するのに有用であることを示した。すなわちピリジニウム様窒素とピロール様窒素に直接結合しているHのNMR化学シフトは0.8と10.8ppmと見積もられた。我々が考慮した様々な端部のN NMRシグナルはピリジン様窒素を除き欠陥でのものとほぼ同じであった。一方、アームチェアー端と欠陥サイトでのピリジン様窒素はその凝集の程度により敏感にN NMRの化学シフトが変化することがわかった。
中島 宏; 中根 佳弘; 増川 史洋; 松田 規宏; 小栗 朋美*; 中野 秀生*; 笹本 宣雄*; 柴田 徳思*; 鈴木 健訓*; 三浦 太一*; et al.
Radiation Protection Dosimetry, 115(1-4), p.564 - 568, 2005/12
被引用回数:8 パーセンタイル:49.16(Environmental Sciences)大強度陽子加速器計画(J-PARC)では、世界最高出力の高エネルギー加速器施設が建設されている。そこで、施設の合理的な遮蔽設計を行うために、J-PARCの遮蔽設計では、簡易計算手法と詳細計算手法を組合せた設計手法が使われている。ここでは、J-PARCの遮蔽設計にかかわる研究の現状について報告する。
安藤 俊就*; 木津 要; 三浦 友史*; 土屋 勝彦; 松川 誠; 玉井 広史; 石田 真一; 小泉 徳潔; 奥野 清
Fusion Engineering and Design, 75-79, p.99 - 103, 2005/11
被引用回数:1 パーセンタイル:10.45(Nuclear Science & Technology)トカマク国内重点化装置のTFコイルはNbAl導体を用い、リアクト・アンド・ワインド法で製作することが検討されている。その製作方法の妥当性を実証するために、実機サイズの導体に0.4%の曲げ歪みを加えて巻き線したコイルサンプルを試作し、試験した結果、その臨界電流値は曲げ歪みを加えていない導体サンプルの臨界電流値とほとんど同じであった。このことは、コイル製作時の曲げ歪みによって臨界電流値が約10%低下するとの予想に反し、実機コイルの製作には朗報である。この現象について解析を行った結果、導体のケーブル部を構成するNbAl線が曲げ加工中に互いに滑ったことが予想された。その考察、解析について議論する。